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STEADY-STATE AND STABILITY ANALYSIS OF EXTERNALLY PRESSURIZED 
GAS -LUBRICATED JOURNAL BEARINGS WITH HERRINGBONE GROOVES 

by David P. Fleming 
Lewis Research Center 

SUMMARY 

A small eccentricity analysis was performed to predict the load capacity and stabil- 
ity of an externally pressurized bearing with herringbone grooves. Numerical results 
were obtained for a range of feeding parameters, pressure ratios, groove lengths, and 
orifice recess volumes for compressibility numbers from 0 to 50. These results were 
obtained from the digital computer program which is presented in an appendix. Results 
showed that the addition of herringbone grooving to an externally pressurized bearing in- 
creases stability, but reduces load capacity at low compressibility numbers. A fully 
grooved bearing is more stable than a partially grooved bearing. Load capacity and 
stability decrease near the speed where the pressure due to the herringbone groove pump- 
ing equals the external supply pressure. Orifice recesses decrease stability, particular- 
ly at high compressibility numbers. 


INTRODUCTION 

Gas-lubricated bearings may be divided into two broad classifications: self-acting 
and externally pressurized. In a self-acting bearing, the film pressure which supports 
the load is developed by the relative motion of the bearing parts and is proportional to the 
fluid viscosity. When there is no motion, the load capacity is zero. In contrast, in an 
externally pressurized bearing, lubricant gas under pressure is supplied from an exter- 
nal source. Thus, this type of bearing can have a substantial load capacity even when it 
is stationary. 

Because of the low viscosity of gases, self-acting gas-lubricated bearings will carry 
a much smaller load than oil-lubricated bearings. For the same reason, they are much 
more susceptible to self-excited instability, commcinly known as fractional frequency 
whirl. A major part of the research in gas-lubricated bearings has been directed toward 



development of bearing configurations that will operate stably. Some of these designs, 
for example, tilting pad bearings, achieve stability at the expense of steady-state load 
capacity. One type of self-acting bearing that has good stability and can also carry a 
higher load than a plain bearing is the herringbone grooved bearing (refs, 1 to 3). 
Externally pressurized bearings also have a higher load capacity than plain self-acting 
bearings and also are fairly stable (ref. 4). 

The principal disadvantage of the externally pressurized bearing is the need for a 
continuous supply of pressurized gas. The herringbone bearing, on the other hand, needs 
no external supply, but has no load capacity at zero speed. The two bearing types could 
be combined; for example, a herringbone grooved rotor could be installed in an external- 
ly pressurized bearing. External pressurization could be used for a startup; upon 
reaching operating speed the external supply could be shut off, and the unit operated as 
a self-acting herringbone bearing. Alternatively, the external supply could be main- 
tained; the inward pumping of the herringbone grooves would reduce the amount of gas 
needed from the external supply . 

Previous analyses have evaluated the load capacity and stability of herringbone 
grooved bearings (refs. 1 and 2) and of externally pressurized bearings (refs. 4 to 6). 
Vohr and Chow determined the load capacity of herringbone grooved bearings in ref- 
erence 1; their analysis was used to evaluate stability in reference 2. Experimental 
stability data for herringbone bearings were obtained in reference 3. The data showed 
that the analysis predicts the onset of instability consistently; however, actual insta- 
bility occurred at somewhat lower speeds than predicted. 

The load capacity of externally pressurized bearings was determined by Lund in 
reference 5. In a later report (ref. 6), Lund calculated the stability of externally pres- 
surized bearings operating at finite eccentricities. Here he included the effect of orifice 
recess volume and attempted to account for there being a finite number of orifices, 
rather than the line source usually assumed. Reference 4 evaluated the stability of an 
unloaded externally pressurized bearing and included the effect of orifice recess volume 
as in reference 6. 

All of these analyses are similar in that they use a small eccentricity perturbation 
and solve for the perturbed pressure using a separation of variables scheme. Thus it is 
easy, in principle, to combine the solutions to find the load capacity and stability of a 
herringbone bearing with external pressurization. 

The objectives of this investigation are to determine analytically the steady-state 
and stability characteristics of externally pressurized herringbone grooved bearings. 
Various combinations of supply pressure, feeding parameter, orifice recess volume, 
and groove length will be explored. 
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ANALYSIS 


The bearing configuration to be analyzed is shown in figure 1. It consists of a double 
row externally pressurized bearing with herringbone grooves. The grooves are shown 
on the rotor, but the analysis is unchanged if they are on the bearing. The herringbone 
grooves may be partial, as shown, or they may extend the full length of the bearing. For 
simplicity of presentation, it will be assumed that the rows of orifices are closer to the 
midplane of the bearing than are the herringbone grooves. This includes the limiting 
case of the orifice rows coinciding with the ends of the grooves. Extension to other cases 
is straightforward. It will be further assumed that the bearing is symmetric about the 
midplane. With this assumption, only half the bearing need be analyzed. Other assump- 
tions are that the number of herringbone grooves is large (ref. 1) and that there are 
enough orifices so that each row may be approximated by a line source. 



Figure L - Externally pressurized herringbone bearing. 


The analysis of reference 1 applies, with the exception that the axial mass flow is no 
longer zero as it was in the herringbone bearing without orifices. Thus, the differen- 
tial equations and boundary conditions must be modified to account for axial mass flow in 
the bearing and flow through the orifice. 

To conveniently obtain solutions for steady whirling, which are needed for the sta- 
bility analysis, a rotating coordinate system is introduced by 
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9* = 9 - u> t 

P 


(1) 


where cu is the frequency of steady circular whirling. Symbols are defined in appen- 

Jr 

dix A. The differential equation to be solved (from ref. 1), is 

| ~ [-*Br sin P + ^zg "^zr> a cos ^ - P Rwh r sin [“^zg + ^ “ “K r ] sin 0 

+ (u) - 0 ) ) — [pah - p(l - a)h 1 sin /3 = 0 (2) 

P gg* L & rj 

This and subsequent expressions were derived for the herringbone grooved section of the 
bearing, but may also be used for the smooth section by setting h^ = h r . The expres- 
sions for the mass flows Jl a . Jl and are in appendix B. The procedure now 

t/r zg zr 

is to approximate the film pressure according to 

p(9*, z) ~ p Q (z) + ep 1 (e*, z) (3) 

Equation (3) and the expressions for the mass flows from appendix B are substituted into 
equation (2). The resulting expression is considered an identity in the eccentricity ratio 
e, and a separate equation Written for each power of e which appears. Powers of e 
higher than 1 are neglected; thus, two equations result. The zero-order equation can be 
written 

£ [“^gO + - a Kro] - 0 < 4 > 


This may be integrated once immediately to yield 

a^zgO + ( 2 " a )^ zr o = ^zO = Constant ( 5 ) 

In terms of the dimensionless pressure Pq = Pg/P a an( * dimensionless coordinate 
£ = z/L, equation (5) for an isothermal bearing becomes 
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The first-order equation is 
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where the Vs are constants given in appendix B. They differ from the constants given 
in reference 1 because dP^/d^ is not constant in a bearing with orifices. The constant 
A is the standard gas bearing compressibility number, and a is a dimensionless 
representation of the whirl frequency. 

Boundary conditions. - At the ends of the bearing 


P 


0 


= 1 


and 
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1 


= 0 


( 8 ) 


At the bearing midplane, by symmetry, 


dp^aPi 

d? a? 


0) 


Pressures Pq and P^ are continuous throughout the bearing film, but there will 
be discontinuities in the derivatives dPg/d? and 9Pj/3£. These are caused at £ = £ , 
by the end of the herringbone groove pattern, and at C = Cj , by the gas flow through the 
orifices. 

Since, in differential equation (6) for Pq, dPg/d£ is the highest order derivative, 
the discontinuity in dPg/dC at £ = £ need not be found explicitly. The discontinuity 
in aP-j/a? at £ = may be determined by noting that Jt^ = otJT ^ + (1 - a )^ zv i is 
continuous at this point. The perturbed mass flows Jt and ^ zr i may be found by 
differentiating the expressions for Jt z and Jt zv (appendix B) with respect to e and 
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then setting e = 0, for example, ^ zgl = 9 ~^ Z g/ 9e l e=0 ° The result for aPj/9? is 


a s 


= ^0^Y7V cose * + *l^ 
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( 10 ) 




g_ 


in which denotes a value of £ infinitesimally greater than £ g . The expression for 
^ , as well as for the other ^'s, is in appendix B. 

At the orifices, conditions are similar to those in an ungrooved externally pressur- 
ized bearing, analyzed in reference 5. One important difference is that, in an externally 
pressurized herringbone bearing, the gas flow through the orifices can be in either direc- 
tion, depending on the supply pressure and the pumping in the herringbone grooves. 
Boundary conditions on pressure are found by balancing the gas flow through the orifices 
with that through the bearing film. The results are 
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The dimensionless mass flow m is given by the usual orifice flow equations, with in- 
herent compensation effects accounted for as in reference 6. 


4R-^TT 


m = Jl r 
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The flow when the bearing is concentric, is obtained by using the concentric values 
P c = P 0c and h f = C in equations (13) and (14). The feeding parameter 
(eqs. (11) and (12)) is defined by 



6^ ViT 
Pa c 3 ifTT?. 


and (in eq. (12)) by 


NV . L _ vL 
jtDLC Pq c I> Pq c d 


(15) 


(16) 


Solution of differential equations . - The constant in equation (6) is determined by the 
gas flow through the orifices. From equation (13), with e = 0, 



Nfp a m 


0 



Combination of equations (15) and (17) with equation (6) gives 


(17) 
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( 18 ) 



m 0 A t 

2P 0*1 


This differential equation contains two unknowns: the pressure P Q and the mass flow 
m Q . The procedure for determining Pq and m Q is (1) assume a value for m^, (2) in- 
tegrate equation (18) from S = 0 to S = Cf, using a forward integration scheme such as 
Runge -Kutta (for ( ^ Cf ? *^^20 = ^ and P 0 =P 0c)’ (3) calculate a new value of nig from 
the value of Pq = Pq(^) just found, and (4) compare with the previous m^; if different, 
repeat steps 2 to 4 until convergence is obtained. 

Differential equation (7) for P ^ may be solved (as in ref. 2) by assuming that P ^ 
is represented by 

P 2 = /?e[G(C)e ie *] (19) 


When this is substituted into equation (7) there results 
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Boundary conditions for Pj in equations (8), (9), (10), and (12) transform into the fol- 
lowing boundary conditions for G 
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The perturbed pressure Pj may now be determined by solving differential equation (20) 
numerically by the method detailed in reference 4. 

Determination of load and stability . - The radial and tangential components of the 
bearing load are found by integrating the film pressure over the bearing area. 


L /'2ir 


H l 


\ -cos 0*1 

p_P \ > R d0* dz 

a l sin 0* J 


Substitution of equations (3) and (19) for P and performance of the 0* integration yields, 
in dimensionless variables, 


rl (Re\ fV 2 


(25) 


The dimensionless forces in equation (25) are defined by 


f r = 


f t = 


ep a LD 


ep a L D 


(26) 


The resultant bearing load W and attitude angle cp may now be calculated. 




1/2 


£ Pa LD 


(27) 


q> = tan 




(28) 


Figure 2 illustrates the relations among these quantities. 

When the bearing is operating stably, the frequency number a is zero. To deter- 
mine the threshold of instability, cr is varied until f^ = 0 (ref. 7). The bearing neutral 
stability condition is then found by equating the centrifugal force, due to the whirling 
bearing mass, to the radial bearing force. 
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Attitude 



Meo; = F 
pn rn 


The subscript n denotes the condition where f t = 0. 
A dimensionless bearing mass may be defined by 


- M P a /C\ 5 
M = — ( — ) 

2Lp 2 ' R ' 


(29) 


(30) 


In terms of previously calculated quantities, M for the neutral stability condition is 
given by 



144 f 


rn 


(31) 


Reference 7 shows that M n is an upper limit of M for stability if the quantity 3f t / da 

is negative at a = a n ; conversely, M n is a lower limit for stability if 3f^/3cr is positive 

at a = or . 
n 


RESULTS AND DISCUSSION 

The analysis of the preceding section has been used to obtain steady-state and sta- 
bility information for a number of externally pressurized herringbone bearing configura- 


10 



tions. Results were obtained using the computer program presented in appendix C. The 
computer program was checked by running cases for a plain externally pressurized bear- 
ing (without grooves) and a herringbone grooved bearing (without orifices). Results for 
the test cases agreed well with those of reference 1 and 4. 

Because of the large number of parameters that may vary in a bearing, the effects 
of all of them were not investigated. Rather, a number of the parameters were fixed. 

The basic bearing chosen for study has a length to diameter ratio of 1 with a single row 
of orifices at the bearing midplane. No inherent compensation effects were included. 

The herringbone groove angle /3 is 30°; the groove width fraction a, 0.5; and the 
groove clearance to land clearance ratio H, 2.1. These values approximate the optima 
found in reference 1 for maximizing the radial load component. Three groove length 
fractions L /L were investigated: 0 (ungrooved), 0.5, and 1 (fully grooved). The feed- 
ing parameter A t was varied from 0 (no orifices) to 4, and the supply pressure ratio 
from 1 to 5. Any other configurations of interest can be easily investigated using the 
computer program in appendix C . 


Steady State Results 

Figure 3 shows the effect of the feeding parameter A^ on load capacity. Figure 3(a) 
is for a pressure ratio P of 1, which means the bearing is actually unpressurized. 

For this case, the load capacity W is greatest when there are no orifices (A^ = 0); W 
decreases with increasing A^. A partially grooved bearing will carry a higher load than 
a fully grooved bearing. The ungrooved bearing is not shown for this case; it is unsuit- 
able for most uses because it is unstable when not loaded. 

Figure 3(b), for a pressure ratio of 2, shows that at low values of A (less than 
about 5) the load capacity now increases with increasing feeding parameter A^.. At 
higher A’s the order is reversed for the grooved bearings. That is, load capacity de- 
creases with increasing A t , as was the case for the unpressurized bearing. The load 
capacity of the partially grooved bearing again exceeds that of the fully grooved bearing. 
For A less than 13 to 33 (depending on feeding parameter) the ungrooved bearing has 
the highest load capacity. At higher A, the load curves for the ungrooved bearings 
level off, while those for the grooved bearings continue to rise. This, of course, is 
because of the increasing self -pressurization by the inward pumping herringbone grooves. 

The load curves for the grooved bearings show an interesting phenomenon in that 
they have a pronounced depression at an intermediate value of A. This depression 
occurs when the pumping of the herringbone grooves raises the pressure on the bearing 
side of the orifices to the pressure that is supplied externally. Near this point, the 
derivative of mass flow with respect to bearing orifice pressure P c becomes very large 
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Pressure 



(a) Feeding parameter, 1. (b) Feeding parameter, 2. (c) Feeding parameter, 4. 

Figure 4. - Effect of pressure ratio on load. Length to diameter ratio, 1; orifice location parameter L f /L, 1; groove angle, 30°; groove width fraction, 0. 5; groove to land clearance ra 




(eq. (14)). That is, a very small change in P c causes a large change in orifice flow. 
Consequently, the bearing is not as well compensated and the stiffness is reduced. 

It should be pointed out that for an actual bearing the loss of load capacity will not 
be as great as predicted by figure 3(b). This is because the difference Am/A(P c /P g ) 
for finite A(P c /P g ) does not approach the infinite value of the derivative 3m/3(P c /P g ). 
Also, the orifice flow equations (14) neglect viscous effects. These become significant 
at low flow rates and reduce the value of 3m/3(P„/P J near P„ = P c . 

Load curves for a pressure ratio of 5 are plotted in figure 3(c) . The trends of fig- 
ures 3(a) and (b) are continued here. Load capacity increases with increasing out to 
the highest compressibility number plotted. No depressions occur in the curves (as in 

fig. 3(b)), because the A value where P„ = P_ is beyond the boundary of the figure. 

c s 

Figures 3(b) and (c) show that the addition of grooves to an externally pressurized 
bearing lowers the load capacity at low compressibility numbers, but at higher A the 
load can be increased. The compressibility number, where the load capacity of the 
grooved bearing first becomes greater than that of the ungrooved bearing, varies with 
the pressure ratio, length of grooves, and feeding parameter. 

The effect of pressure ratio on load is shown in figure 4. Load capacity generally 
increases with pressure ratio. The load capacity for P a = 2 is little different than for 
P g = 1 (unpressurized), particularly at higher compressibility numbers. Near the value 
of A where P = P (depression in load curve), the load capacity for P = 2 can drop 
below that for P g = 1. Increasing the pressure ratio to 5 results in a relatively large 
increase in load capacity, particularly at low compressibility numbers and large feeding 
parameters . 

Attitude angles are plotted in figure 5 for a pressure ratio of 2. Except at quite low 
compressibility numbers (A < 5), attitude angles are smaller for smaller values of the 
feeding parameter A^.. The grooved bearings (figs. 5(b) and (c)) exhibit behavior gener- 
ally similar to an ungrooved externally pressurized bearing, with two exceptions. At 
zero speed (A = 0) the grooved bearings have a small attitude angle, which may be posi- 
tive or negative, depending on the length of the grooves. In addition, when the pumping 

of the herringbone grooves increases the pressure P to near the supply pressure P , 

c s 

the attitude angle rises rapidly. This corresponds to the drop in load capacity mentioned 
earlier. A comparison of figures 5(a) to (c) shows that attitude angle decreases with in- 
creasing groove length, except at high compressibility numbers (A > 35), where the bear- 
ing with half-length grooves has a lower attitude angle than either the fully grooved or 
ungrooved bearing. 

Because the herringbone grooves act as a pump, the gas pressure in the bearing in- 
creases with increasing compressibility number. Figure 6 shows the pressure at the 
center of the bearing, Pq c , for three values of the feeding parameter A t , at an external 
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Figure 5. - Effect of feeding parameter on attitude angle. Pressure ratio P< ■% length to diameter ratio, 1; orifice location Figure 6. - Pressure at center of bearing. Pressure ratio, 

parameter L f /L, 1; groove angle, 30°; groove width fraction, 0.5; groove ft land clearance ratio, 2.1. 2; length to diameter ratio, 1; orifice location parameter 

Lf/L, 1; groove angle, 30°; groove width fraction, 0. 5; 
groove to land clearance ratio, 2. 1. 



supply pressure ratio of 2. As expected, larger values of A^. decrease the rate of 
change of Pq c with A. This is because a larger A^. means the bearing is more open 
to the external supply. The fully grooved bearing’s pressure changes more rapidly than 
does the partially grooved, since pumping is proportional to groove length. Pressure 
does not change with compressibility number in the ungrooved bearings; to avoid confu- 
sion, these curves have not been extended to A=0. The solid lines, for A t = 0, indi- 
cate the pressure in a bearing with no orifices; they are, of course, valid for any supply 
pressure. 

The gas mass flow through the bearing appears in figure 7 for the same bearing con- 
ditions as in figure 6. There are additional curves for A^. = 1. For A^ = 0, the mass 
flow is always 0. A negative value of mQ indicates the gas flow is from the pressurized 



Figure 7. - Gas flowthrough bearing. Pressure ratio P s , 
2; length to diameter ratio, 1; orifice location parameter 
Lf/L, 1; groove angle, 30°; groove width fraction, 0. 5; 
groove to land clearance ratio, 2. 1. 
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supply to the bearing, whereas positive m Q indicates the bearing is pumping gas into 
the pressurized supply. 

The mass flow curves are arranged in the same order as the bearing pressure curves 
of figure 6, Note that the magnitude of the dimensionless flow mQ decreases as the 
feeding parameter increases. This is due to the way mass flow is made dimensionless. 

At zero speed (A = 0), gas consumption increases with increasing length of grooves. 

This is because a grooved portion of a bearing has a larger flow area, and thus offers 
less resistance than an ungrooved portion. 


Stability Results 

In order to keep the figures presented to a reasonable number, stability information 
will not be given for all combinations of groove length, pressure ratio, feeding param- 
eter, and orifice recess volume which were investigated. Instead, the basic bearing 
mentioned at the beginning of this section will be further defined, and results presented 
for variations of each of the preceding four parameters from their basic values. These 
basic values are L /L = 1, P = 2, A, = 2, and v = 0. 

o S C 

Figure 8 shows the variation of stability with compressibility number, with groove 
length as a parameter. Stability, as measured by the dimensionless mass M, generally 
decreases with increasing A, and increases with groove length. Above a compressibility 
number of 15, the stability of the fully grooved bearing increases sharply and becomes 
much greater than that of either the half-grooved or ungrooved bearing. 

The stability curves of the fully grooved and half-grooved bearings have distinct de- 
pressions near A = 14 and A = 27, respectively. These depressions correspond to the 
depressions in the load curves of figures 3(b) and 4(b). As was discussed regarding the 
load capacity curves, the drop in stability in an actual bearing whirling with a finite 
eccentricity would probably not be as drastic as figure 8 predicts. With an actual bear- 
ing, problems with this low stability region can probably be avoided by passing through 
the region rapidly, either by accelerating the bearing rotor or changing the supply pres- 
sure. 

Near a compressibility number of 38, the stability curve for the half-grooved bear- 
ing becomes very steep, in effect imposing an upper speed limit even for very small 
values of M. This limit generally occurred between compressibility numbers of 20 and 
40 in the half-grooved bearings. A similar limit was not observed for fully grooved or 
ungrooved bearings within the range of compressibility numbers investigated. 

The effect of the feeding parameter A^. on stability is shown by figure 9. At low 
compressibility numbers (A < 5), higher feeding parameters give greater stability. 

Near A = 14, where the herringbone pumping pressure becomes equal to the supply 
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Figure 8. - Effect of groove length on stability. Pressure ratio P s , 2; feed- Figure 9. - Effect of feeding parameter on stability. Pressure ratio P s , 

ing parameter 2- orifice recess volume ratio, 0; length to diameter ratio, 1; groove length ratio LJL, 1; orifice recess volume ratio, 0 ; length bn 

orifice location parameter L f /L; 1; groove angle, 30° : groove width frac- ameter ratio, 1; orifice location parameter L f /L, 1; groove angle, 30 ; 

tion, 0. 5 : groove to land clearance ratio, 2. 1. groove width fraction, 0. 5; groove to land clearance ratio, 2. 1. 



pressure, the order is reversed, with the no-orifice bearing (A^. = 0) most stable. At 
high A (>20), there is no clear trend. The greatest stability is offered by A^ = 2, and 
the least by A^ = 4. 



Figure 10. - Effect of pressure ratio on stability. Groove length ra- 
tio, L/L, 1; feeding parameter, A t , 2; length to diameter ratio, 1; 
orifice recess volume ratio, 0; orifice location parameter Lf/L, 1; 
groove angle, 30°; groove width fraction, 0.5; groove to land 
clearance ratio, 2. 1. 


Figure 10 shows the effect on stability of a variation in supply pressure ratio. In- 
cluded in this figure is the curve for A^. = 0; the remainder of the curves are for A^ = 2. 

For low compressibility numbers (A < 10), P = 5 provides the greatest stability; the 

s 

curve is beyond the maximum ordinate of the figure. Stability decreases with decreasing 

pressure ratio. The plain herringbone bearing (A^. = 0) is between P g = 1 and P g = 2 

(P_ = 1 denotes a bearing whose supply lines are open to the atmosphere). For compres- 

sibility numbers between 10 and 20 there is no clear optimum, as the curves rise and 

fall. At high A, pressure ratios of 1 and 2 are more stable than P = 5 or the no- 

s 

orifice bearing. The figure shows that the addition of orifices to a herringbone bearing, 
without pressurization, lowers the stability at low compressibility numbers, but can in- 
crease the stability at high A . 

The decrease in stability due to a small orifice recess is shown in figure 11. For 
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Orifice 

recess 

volume 

ratio, 

v 

0 

2 


Figure 11. - Effect of recess volume on stability. Pressure ratio, 2; groove 
length ratio, L„/L, 1; feeding parameter, 2; length to diameter ratio, 1; 
orifice location parameter L f /L, 1 ; groove angle, 30°; groove width frac- 
tion, 0.5; groove to land clearance ratio, 2.1. 


compressibility numbers less than 14 there is only a small loss of stability. At higher 
A, however, the stability for a recess volume ratio of 0.2 continuously decreases, while 
the stability for v = 0 increases from A = 14 to a maximum at A = 40 before decreas- 
ing. This behavior at higher compressibility numbers is typical of the effect of recess 
volume on the bearing configurations studied. 

Multibranch curves . - For all cases of finite recess volume in grooved bearings, 
there was more than one neutral stability condition found at the higher compressibility 
numbers. That is, for a given A, there was more than one whirl frequency which yielded 
f^. = 0. Figure 12 illustrates this for a feeding parameter of 2, pressure ratio of 5, and 
recess volume ratio of 0.2. The controlling curve will be the lowest; this will give the 
maximum value of M at which the bearing will be stable. Therefore, in using the anal- 
ysis to determine stability, one must use some caution to be certain the smallest value of 
M n has been found. 
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Figure 12. - Multibranch curves. Pressure ratio, 2; feeding parameter, 2; ori- 
fice recess volume ratio, 0. 2; length to diameter ratio, 1; orifice location 
parameter L f /L, 1; groove angle, 30°; groove width fraction, 0.5; groove to 
land clearance ratio, 2. 1. 


SUMMARY OF RESULTS 


A small eccentricity analysis was performed to determine the steady-state and sta- 
bility characteristics of externally pressurized bearings with herringbone grooves. 
Compressibility numbers from 0 to 50 were investigated. The following results were 
obtained for a bearing with a single row of orifices and a length to diameter ratio of 1: 

1. The addition of herringbone grooves to an externally pressurized bearing increases 
the stability of the bearing. Grooving reduces load capacity at low compressibility num- 
bers but increases load capacity at high compressibility numbers. 

2 . The fully grooved bearing is generally more stable than the partially grooved 
bearing. This is especially true at high compressibility numbers. However, the par- 
tially grooved bearing has a higher load capacity. 

3. At low compressibility numbers, stability increases with increasing supply pres- 
sure and feeding parameter. At high compressibility numbers, there is no clear rela- 
tion of stability with feeding parameter and supply pressure. 
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4. Load capacity and stability decrease, and the attitude angle increases near the 
speed where the pressure due to the herringbone groove pumping equals the external 
supply pressure. This at least partly due to assumptions in the analysis. 

5. Orifice recesses decrease stability. The effect is marked at high compressibility 
numbers. When there are orifice recesses, more than one neutral stability condition 
can exist at high compressibility numbers. The controlling condition is that which gives 
the lowest dimensionless mass for neutral stability. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 30, 1970, 

129-03. 
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APPENDIX A 


SYMBOLS 


a 

orifice radius 

M 

C 

ridge clearance at zero eccen- 



tricity 

m 2 

Sf 

coefficient (see appendix B) 


D 

bearing diameter 

m 

d 

orifice recess diameter 

N 

e 

journal eccentricity 

P 

F 

bearing load component 

P 

f 

dimensionless load component, 



F/ep a LD 

Pa 

G 

dimensionless complex function 

R 


of ? 

3t 

H 

ratio of ridge clearance to groove 

(le 


clearance when bearing is con- 

T 


centric, (h g /h r ) Q 

t 

h g 

local film thickness over groove, 


C(H + e cos 0*) 

h r 

local film thickness over ridge, 

u p 


C(1 + e cos 6*) 

V 

9m 

imaginary part of expression 

V 

i 



k 

specific heat ratio 

W 

L 

bearing length 

W 

h 

length of bearing outboard of 

z 

orifices 



total axial length of grooves 

Zf 

M 

rotor mass per bearing 


Jf 

lubricant flow per unit length 

z g 


dimensionless rotor mass, 
(MP a /2L M 2 )(C/R) 5 

dimensionless rotor mass, 
MCa> 2 /p 0 LD 

d 

dimensionless lubrica nt flow rate, 
4^ Z R yfdt T Vl + 5 2 / Na 2 p a 

number of orifices per bearing 

dimensionless pressure, p/p„ 

ci 

pressure 

atmospheric pressure 
bearing radius 
gas constant 
real part of expression 
absolute temperature 


surface speed of grooved member 

surface speed of smooth member 

orifice recess volume 

orifice recess volume ratio, 
NV/jtDLC 

total bearing load 

dimensionless load, W/ep LD 

ci 

axial coordinate measured from 
end of bearing 

distance from end of bearing to 
first row of orifices 

L /2 

W 
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a ratio of groove width to width of 
groove -ridge pair 

orifice discharge coefficient 

/3 groove angle (fig. 1) 

5 inherent compensation factor, 

a 2 /dC 

e eccentricity ratio, e/C 

£ dimensionless axial coordinate, 

z/L 

£ f z f /L 

V L 

0 angular coordinate 

9 * rotating angular coordinate , 

A bearing compressibility number, 

6pa)R 2 /p a C 2 

A A(L/R) if smooth member is 

s 

rotating; -A(L/R) if grooved 
member is rotating 

A t feeding parameter, 

6julN 2 V^T / P a C 3 Vl + 5 2 

/i lubricant dynamic viscosity 

p local lubricant density 

2 2 

a frequency number, 12pc0pR /l^C 

cp attitude angle 


* 0 -(V 2P 0c)[ 3m / S(P c/ P s>]| £=0 

NV/nD 2 CP 0e 
o> rotational speed 

whirl frequency 
Subscripts: 

c condition immediately downstream 

of orifice 

g groove region 

n condition at which f^. = 0 

r radial; ridge region 

s condition upstream of orifice 

t tangential 

z axial direction 

9 circumferential direction 

0 zero eccentricity 

1 perturbed quantity 

Superscripts: 

+ value of coordinate infinitesimally 

greater than base value 

value of coordinate infinitesimally 
less than base value 
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APPENDIX B 


EXPRESSIONS FOR BEARING MASS FLOW AND 
DIFFERENTIAL EQUATION COEFFICIENTS 

Mass Flow Expressions Used in Equation (2) (from ref. 1) 


ji - 
Zg 


12m ahl + (1 - a)h 2 


((1 - a) (hi - h 2 ) sin 0 cos 0 - ^2- 
L v g r/ R 90 


h 2 +(l-a)(h 2 -h 2 ) sin 2 0 


- 6 M (1 - ot)(U - U )(h - h ) sin 0 cos 0 

9Z P & 5 1 


~^zr " 


12 ^ + (1 - a) h| 


j-a (h 2 - h 2 ) sin 0 cos 0 i- 
I v g r/ R 90 


h g ' a ( h g “ h r) sin2 ? 


+ 6|LiO!(U p - U g )(h g - h r ) sin 0 cos 0 


h' 


0r 


12m ah 2 + ( 1 - a)hl L L 


h g ' «( h g " h r) cos2 P 


R 90 


a(hg - h 2 ) sin 0 cos 0_3E - 6 M a-(U p - U g )(h g - h p ) sin 2 pj + ^ (U p + U g ) 


Differential Equation Coefficients 


*1 = 


H 2 + a( 1 - o:)(h 3 - if sin 2 0 R 


a + (1 - a)H v 
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a( 1 - a)(H° - 1>(H - 1) sin 0 cos 0 A c 

« p = s 

H 3 + a(l - a)(H 3 - l) sin 2 0 

<€ = 2a (1 - Q!)(h 3 - l) sin 0 cos 0 L 
H 3 + a(l - a)(H 3 - l) 2 sin 2 0 R 


»4a = ^3 


’4b 


= (gH + 1 - a)[a + (1 - a)H 3 ] /L\ 2 




H 3 +a(l - a)(H 3 - 

n 2 9 \R 
1) sin (3 

a(l - a)(H 3 - 1)(H - 

l)A g sin 2 j3 

H 3 + a(l - a)(H 3 - 

\ 2 2 ^ 
lj snY 0 

. H 3 + a(l - a)(H 3 - if 

2 0 
cos* 0/LY 


H 3 + a(l - a)(H 3 - l) sin 2 0 


^R/ 


'6a 


3a(l - a)(H 3 - 

1)(H - 1) sin 0 cos 0 [a(H 2 - H - 1) - (1 - a)H 2 (H 2 + H - 1) 


2 

H 3 + a(l - a)(H 3 - l) sin 2 0 

[a + (1 - a)H 3 ] 


'6b 


a + (1 - a)H 3 ] 


H 3 + a(l - a)(H 3 - l) sin 2 0 


«6o = 


3a(l - a)(H - l) 2 H 2 A o sin 2 0 

S 


H 3 + q;(1 - q!)(h 3 - l) sin 2 0 [a + (1 - a)H 3 ] 


L 

R 




3a(l - a)(H - 1)(H 2 - l)A g sin j3 cos 0 H 2 (H 2 + 1) + a(l - a)(H 3 - l) sin 2 (3 


2 

H 3 +a(l - a)(H 3 - l) sin 2 0 
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5=0 1 



_ 3aH 2 (H - l){l + (1 - g)(H 3 - 1)[2 + (1 - g)(H 3 - 1)] sin 2 i 

m T 2 l 2 

H 3 +a(l - a)(H 3 - l) sin 2 f3 

In the ungrooved portion of the bearing, these coefficients become 


<fp = ° 

Sf 3 = ° 



= o 



APPENDIX C 


COMPUTER PROGRAM FOR ANALYSIS OF EXTERNALLY PRESSURIZED 
HERRINGBONE GROOVED GAS-LUBRICATED JOURNAL BEARING 

The program for determining steady-state and stability characteristics of the exter- 
nally pressurized herringbone bearing is written in FORTRAN IV, version 13, for use on 
the IBM 7094-11 digital computer at Lewis Research Center. Minor modifications may be 
necessary to allow the program to be used with other computing systems. 


Program Input 


Two variables are set within the program: k, the specific heat ratio of the lubricant 
gas, and RKEP, an accuracy parameter. As the program operates, the step size in the 
Runge-Kutta differential equation solver is successively halved until the magnitude of the 
bearing load W varies by less than RKEP times W. 

The remainder of the information needed by the program is read from.punched cards. 
A triad of cards shows each configuration to be analyzed. The first gives geometric and 
operating parameters; the second contains the array of bearing numbers A to be used 
(up to 19 values); and the third gives initial estimates of the frequency number a, one 
for each A value. Specific formats follow. 

Geometric configuration card : Format (13F6. 0). - Succeeding eight- column fields 
contain, in real format (with decimal point): /3, a, H, L/D, L g /L, SGN, Lj/L, A t , P g , 
v, a d , and 1 +5 2 . SGN is +1.0 if the smooth member of the bearing is rotating, and 
-1.0 if the grooved member is rotating. 

Bearing number array card: Format (14, 19F4.1). - First four -column field: the 
number of A’s in the array, integer format, right adjusted. Succeeding four-column 
fields: the values of A for which calculations are desired, real format (with decimal 
point). Positive A means the program will search for a root of f t = 0 such that 
9fj./3cr<0. The reverse is true if A is negative. 

Frequency number estimate card : Format (4X, 19F4.0). - First four- column field 
is not read. Succeeding four- column fields contain, in real format (with decimal point), 
the initial estimate of frequency number a, one value for each A. 

Any number of geometric configurations may be examined in one run of the program. 
Sample program input on a FORTRAN coding sheet appears as figure 13. 

The first card is a geometric configuration card which indicates that j3 = 30°, 


a = 0. 5, H = 2. 1, L/D = 1, L g /L = 1, SGN = -1, h f /L = 1, A t = 2, P g = 2, v = 0, a d = 0, 
and 1 + 5 2 = 1. The second card, the A array, indicates there are 13 A values ranging 
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STATEMENT 

NUMBER 

1 

FORTRAN STATEMENT 

IDENTIFICATION 

12 3 4 5 


7 B 9 10 II 12 

IS 14 15 16 17 16 

19 2 0 21 22 2 3 24 

ZS 26'2? 28 29 30 

31 32 33 34 3S 36 

37 38 39 40 41 42 

45 44 45 46 47 4B 






■EE 


. 5 

■MHW 

1 . 

1 . 

-1 . 

1 . 


mmm 

o. 

0. 

mm 


13 


mm 


2 . 4 . 

8 . 1 



0. 2 4. 

mm 

0 . 50 . 


■ 




spra 



5 . 

MMfSfl 

nsMfl 

0. 5 . 

5 . 

5 . 5 . 










mm 

MM 



MH 










mm 

mm 










‘ ' 

|fgg|g 


MM 

■ 


N N 




ippm 





H 








■Ml 

■PH 




, , , , , 






m 



■ 

MPP| 













jjjp 










mm 



ppm 

■PPP^ 










N 



ppm 

PPM 










1 

1 mh 







mm 

mm 












n 

m i 













M 

M 

mm 








■ 

pm 














PE 


i 








pm 














pp 

■ 









MM 




pm 







■W 



1 

■ 




M | 







M 























■MM 











ppp 









■ 





i p 




















a 



19 20 21 22 23 24 

25 26 27 28 29 30 

31 52 33 34 35 36 



49 50 51 52 53 54 

55 56 57 58 59 60 



masm 


Figure 13. - Sample program input. 


from 0 to 50. The third card is the a estimate card, showing initial estimates of a, 
which correspond to the A array. Any number of these sets of three cards may be 
used. 

When the orifice discharge coefficient is specified as zero on the geometric config- 
uration card, the actual value used is computed as a function of the pressure ratio across 
the orifice, using the information in reference 2. If a nonzero a< d is specified, that 
value is used. 


Program Output 

Output consists of two sheets for each configuration: A working sheet and a calcu- 
lation summary sheet. 

The working sheet shows the geometric and operating parameters, trial values for 
determination of the pressure downstream of the orifices, final value of orifice down- 
stream pressure, i^q, andRKEP. For each A in the array and for each trial value 
of a there appear A, f r , f^, W/(ep a LD), cp, Mg, final Runge-Kutta step size, M, 
Wp/w, a, and a coefficient applied to if/ Q in equation (24). This coefficient, in refer- 
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ence 4, was taken as 1 when a was 0, and 0 otherwise. In the present work it is always 
1. A sample working sheet appears as figure 14. The data appearing thereon were gen- 
erated from the values shown in the sample input. 


EXTfc RNALLY PRESSLRIZFD HERRINGBONE BEARING — GROOVED MEMBER ROTATING 


BETA 

ALFA 

H 

L/D 

LIGJ/L 

L (0)/L 

LAMBDA T 

PS/PA 

VC 

CRF.COEF 

l+0» «2 

K=CP/CV 

30. CO 

0.500 

2. IOC 

1.000 

1.000 

1.000 

2.000 

2.000 

0 

VARIABLE 

1.C00 

1.400 


SCLUTICN 

FOR 1 

POC LAMBDA* C 













1 

PC 

1.4142136 PCI 

1.3121654 

MO -0.8383 


MAXMO 

0 


UNCHOKED 

OR IF 

CQEF 0.661 




2 

PC 

1.3631895 PCI 

1.3026086 

MO -0.8687 


MAXMO 

0 


UNCHOKED 

CRIF 

COEF 0.67U 




3 

PC 

1.3026086 PCI 

1.2983005 

MO -0.9015 


MAXMO 

0 


UNCHOKED 

CRIF 

CQEF 0.682 




A 

PC 

1.2983005 PCI 

1.2982814 

MO -0.9037 


MAXMO 

0 


UNCHOKED 

CRIF 

COEF 0.682 



PC 

1 . 

298 

PC/PS 0.649 

PSO-O. 783 

PS1 0 


RKDX 

REDUCED 

UNTIL ERROR 

.LE. 

0.500E-02 



LAMBOA 

FR/EPALO FT/EPALO 

W/EPALD 

PHI MCW2/PAL0 


RKDX 

MPC5/2LR5MU2 WHIRL RATIO V 

IE NR SIGMA 

PSC COEF 

0 


0 

.282 -0.701E-C2 

0.2819 

-1.425 

0 



0.1250 

0 


10CC.0 

C 

l.COOO 

0 


0 

.291 0.775E-C1 

C. 3012 

14.90 

c 



0.1250 

41.916 


1C00.0 

-l .CCCO 

1.C000 

0 


0 

.282 0.919E-C4 

0.2818 

0. 187E-01 

0 



0.1250 

5091.0 


10C0.0 

-0.8 300E-0 1 

l.COOO 


SCLUTICN 

FOR POC 

LAMECA= 

0.500 










l 

PC 1.4142136 PCI 

1.3337132 

MO 

-0.8383 

MAXMO 

0.9587E-01 

UNCHOKED 

CRIF 

COEF 0.661 



2 

PC 1.3739634 PCI 

1.3285685 

MO 

-0.8625 

MAXMO 

0.9587E-01 

UNCHOKEO 

ORIF 

COEF 0.668 



3 

PC 1.3285685 PCI 

1.3264548 

MO 

-0.8878 

MAXMO 

0.9587 E-0 1 

UNCHOKED 

OR I F 

COEF 0.677 



4 

PC 1.3264548 PCI 

1.3264503 

MO 

-0.8890 

MAXMO 

0.9587E-0 l 

UNCHOKEO 

CRIF 

COEF 0.677 



PC 1. 

326 PC/PS 

0.663 

PSO-O. 807 


P$1 0 

RKDX REDUCED 

UNTIL ERROR 

.LE. 

0 • EOOE-02 



LAMBDA 

FR/EPALO 

FT/EPALD 

W/EPALD 

PHI MCW2/PALD 

RKDX 

MPC5/2LR5MU2 WHIRL RATIO VI 

te \R SIGMA 

PSC COEF 

0.5CC 

0.313 

0.27 1E-C1 

C. 3137 

4. 

959 

0 

0.1250 

0 


0 

0 

1 .0000 

0.5CC 

0.311 

0.106E-C1 

C • 31 14 

1. 

956 

7.7793 

0.1250 

1120.2 


0.2000 

C.2C00 

1.0000 

0.500 

0.311 

-0.594E-C2 

C. 3106 

-1. 

096 

1.9410 

0.1250 

279.51 


0.4C00 

0.4CC0 

l.COOO 

0.500 

0.311 

0.58 1E-C5 

C • 3107 

0. 

107E-02 

2.8828 

0.1250 

415. L2 


0.3283 

C.3283 

l.COOO 


SCLUTICN 

FOR POC LAMECA* 1 

L.COO 









I 

PC 1.4142136 PCI 

1.3576819 

MO -0.8383 

MAXMO 

0.1917 

UNCHOKEO 

CRIF 

COEF 0.661 



2 

PC 1.3859477 PCI 

1.3554464 

MO -0.8554 

MAXMO 

0.1917 

UNCHOKEO 

CRIF 

COEF 0.666 



3 

PC 1.3554464 PCI 

1.3545969 

MO -0.8731 

MAXMO 

0.1917 

UNCHOKEO 

ORIF 

COEF 0.672 



4 

PC 1.3545969 PCI 

1.3545963 

MO -0.8735 

MAXMO 

0.1917 

UNCHOKED 

ORIF 

COEF 0.672 



PC 1. 

355 PC/PS 0.677 

PSO-O. 830 

PSL 0 

RKDX REDUCED 

UNTIL ERROR 

.LE. 

0. 500E-02 



LAMBDA 

FR/EPALD FT/EPALD 

W/EPALD 

PHI MCW2/PALD 

RKDX 

MPC5/2LR5MU2 WHIRL RATIO VIE 

NR SIGMA 

PSC COEF 

1.000 

0.347 0.575E-C1 

0.3513 

9.414 

0 

0.1250 

0 


0 

0 

1 .COOO 

l.COO 

0.339 -0 .22 CE-Cl 

C • 3394 

-3.718 

1.3547 

0.1250 

48.769 


0.5000 

l.CCCO 

l.COOO 

l.COO 

0.34C 0.101E-C1 

C • 3398 

1.708 

3.7743 

0.1250 

135.87 


0.3000 

C .6C00 

l.COOO 

1.000 

0.339 0.252E-C4 

0.3391 

0.426E-02 

2.5726 

0.1250 

92.613 


0.3630 

0 . 7261 

l.COOO 


SOLUTICN FOR 
1 PC 

POC LAMBDA* 

1.4142136 PCI 

2. COO 
1.4108198 

MO 

-0.8383 

MAXMO 

0.3835 

UNCHOKEO 

CRIF 

CCEF 

0.661 

2 PC 

1.4125167 

PC 1 

1.4108133 

MO 

-0.8394 

MAXMO 

0.3835 

UNCHOKEO 

CRIF 

COEF 

0.661 

3 PC 

1.4108133 

PCI 

1.4108112 

MO 

-0.8404 

MAXMO 

0.3835 

UNCHOKED 

ORIF 

COEF 

0.662 

PC 1.411 

PC/PS 0.7C5 


PSO-O. 874 


PS1 o 

RKDX REDUCED 

UNTIL ERROR 

.le. 

0. 500E-02 


Figure 14. - Example of working sheet. 


The calculation summary sheet again shows the geometric and operating parameters, 
and, for each A of the A array, W/(ep a LD), tp, c o /w, <r n , M 2n , M n , P Qc , (9f t /9ff) n , 
m Q , oi^ (actual value used), and the ipQ coefficient. Figure 15 shows an example of a 
summary sheet. Again, the data correspond to those of the sample input. Program 
execution time for the sample input was 20 seconds. 
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EXTERNALLY PRESSLR IZEO HERRINGBONE BEARING — GROOVED MEMBER ROTATING 


L/D 

LO/L 

LAMBDA T 

1+DELTA**2 

K=CP/CV 

ORIF COE F REC • VOL. VC 

1.000 

i.ooo 

2.000 

1.000 

1.400 

VARIABLE 0 

PS/PA 

BETA 

ALFA 

H 

LG/L 

RKEP 

2. COO 

30.00 

0.500 

2.100 

1.000 

0.500E-02 


BRG NR 
LAMBDA 

BRG LOAD 
W/EPALO 

ATT ANGLE 
PHI 

WHIRL RATIO 
WP/OMEGA 

VIB NR 
SIGMA 

STABILI TY 
MCW2/PALD 

PARAMETERS 

MPC5/2LR5MU2 

PC/PA 

DIFT ) 
DISIGMA) 

MASS 
FLOW PC 

ORIFICE 
COEF AO 

PSO 

COEF 

0 

0.282 

-1.425 

1000. -0.830E- 

01 0 

0. 589 E 04 

1.298 

-0.84 E— 0 1 

-0.904 

C.682 

1 .000 

0.5C0 

0.314 

4.959 

0.328 

0.328 

2.883 

415.1 

1.326 

-0.83E-01 

-0.889 

C.677 

1.000 

1.000 

0.351 

9.414 

0.363 

0.726 

2.573 

92.61 

1.355 

-0.80E-01 

-0.874 

C.672 

1.000 

2.CC0 

0.434 

14. 10 

0.371 

1..484 

2.863 

25.77 

1.411 

-0.76E-01 

-0.840 

C. 66 2 

1.000 

4 .000 

0.592 

15.29 

0.355 

2.840 

3.948 

8.884 

1.523 

— 0.66E-0 1 

-0.765 

C.642 

1 .000 

8.CC0 

0.815 

11.55 

0.313 

5.005 

6.796 

3.823 

1.743 

-0.56E-01 

-0.572 

C.60E 

1.000 

12.00 

0.920 

13.2 3 

0.340 

8.162 

5.828 

1.457 

1.945 

-0.47E-01 

-0.271 

C. 584 

1 .000 

13. CO 

0.896 

20.24 

0.397 

10.32 

3.318 

0.707 

1.984 

-0.57E-01 

-0.146 

C. 5 8 1 

1.000 

16. CO 

1.022 

9.870 

0.306 

9.778 

8. 340 

1. 173 

2.079 

-0.37E-0 1 

0.329 

C. 586 

1 .000 

20. CO 

1.138 

4.334 

0.183 

7.303 

30.48 

2.743 

2.294 

—0 .18 6—0 1 

0.655 

C.60E 

1.000 

24.00 

1.238 

2.377 
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Figure 15, -Calculation summary sheet. 
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UBFTC EPHB 


C SMALL EOCENE* It I FY ST Ab I L 1 FY ANALYSIS F UK HERR I NS BuNt BEARING w 1 1 H OR [FILLS 
DIMENS ION SLAM(20)»SW(20) * iPH 1(20) , SW3 ( 20 ) , S CrM ( 20 ) , SS I G ( 20 ) 

DIMENSION ACM ( 20) ,SD( 2 0) , SPC (20) ,SMU (20) ,SaU (20 ) ,S PS ( 2D ) 

COMMON /CY4/ DP, DC, DS, OOP 

COMMON /CFM/ PCK ,K , SKI ,SK2 , VARA ,A0 

common /cm/ CC,UC,C 

COMMON /CF F/ HLM,CP ,C 2 A , C2B , I CO A ,IC4A , I C4 W , 1 C SAW , I Co A , I C6W , I Cs 6 W , 

1 C5A, C 7A , C 83 » HLM2 » C S5 

COMMON /CHS/ LAMBDA ,PSO,PSL ,C MA ,CMC , W , PHI , wR, RKEP, Q,L, PC, 0 48,040, 

L C6B,C6C , FUST, SI Z,CPS 
COMMON /CHR/ CG,OD,ZG,ZO 
REAL L,LG,L0,LF,!M0,C(2) ,K,LS, LaMBDA 
COMPLEX 1 C3A, IC4A, IC4W, ICS4W, IC6A, LCoW,ICS6rt 
COMPLEX Y(2),A(2),B(2),XF,CDX,CG,0G 
LOGICAL VARA »F IRSF 
EXTERNAL FI 

ALPHA(P) = .84276 + (. 06461 + (-.76195 f . 4 34 if> *P ) *P ) *P 

EPS = l.E-4 
EPSW = l.E-3 
K - 1.4 

SKI = SQRT (2. *K/(K+l. ))* <2. /(K+l. ))**(!./ (K-l.) ) 

SK2 = SOR T ( 2 . * K/ ( K- 1 . ) i 
PCR = (2./(<+L.))**U/(K-l.i) 

RKEP = .006 

3 READ (5,5) BE TAD, ALF A , H , CL , LG , SON , LO , L T , PS , VC , AO, D l 

5 FJRMAF (13F6.0) 

C BETAD = GROOVE ANGLE, OEGREES. ALFA = GKUUVt WIDTH FRACTION. H = GROJVF 
C CL EAR AN C £ /L AN 0 CLEARANCE, CL = L/D. LG = GROOVE LENGTH FRACTION uG/l. 

C SGN = +1.0 FOR SMOOTH MEMBER ROTATING, -1.0 FOR GROOVED MEMBER ROTATING 
C LO =FR AC T I ON OF LENGTH OUTBOARD OF ORIFICES, LT- L AM BOA( T ) , PS = 3 RESSJRE RATIO 
C VC=ORIFICE RECESS' VOLUME RATIO, AO=QRIFICE COEFFICIENT, 01 = UUELT A**2 
C LG= l DENOTES FULLY GROOVED BEARING, L0=1 DENOTES SlNGl.c ROW OF ORIFICES 
IF (LG.EO .0. ) SGN = 0. 

VARA = .FALSE. 

IF (AO.Ej.O.) VARA = .TRUE. 

READ (5,1) IS, (SLAM(I), I = l, IS) 

I FORMAT (14, 19F4. 1) 

READ ( 5,2 ) ( SSIG( I ) , I = i , IS) 

FORMAT (4X, 19=4.0) 

FORMAT ( 1H1/8X»65HSMALL ECCENTRICITY ANALYSIS FOR HERRINGBONE BEAR 
LING W I TH ORIFICES ) 

PCA = P S* PCR 

BETA = SIGN ( BETAD/57. 2957 /96 , SGN) 

L = 2.*CL 
S 1= S INI 3 ETA) 

C0= CUS(BETA) 

ALFI= I .-ALFA 
H2= H* H 
H 3= H 2*H 

0N= ALFA+ALF I#H3 
H32= (H3-l.)**2 
AA= AL FA* ALF 1 
AS= AA *S I * SI 
ASC= AA*SI*CO 
02= H 3+AS*H 32 
H3l= ( H3- 1.)* (H-l. ) 

Ci= U2/(DN*L ) 

WRITE (6, 707) 
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70/ FORMAT { iHl/1 5X»45HEX TfcRNALLY PRESSURIZED HEKR ING8 0NL BEARING -- ) 
IF I SGN) 9, 9G S 91 
91 WRITE (6, 808) 

808 F JRMA T ( 1H + 60X, 22HSMU0H MEMBER RUT ATI NG ) 

GO TO 902 

90 WRITE (6,890) 

890 FORMAT ( 1H + 60X, 1 7rlUN3ROO VEO BEARING ! 

GU TU 902 

9 WRITE (6, 909) 

909 FORMAT ( 1H+ BOX, 23H3RUUVE0 MEMBER ROTATING ) 

90 2 WR I TE ( t>, 6 ) 

6 FORMAT ( 6 HL 3ETA 7X.3HALFA 8X,LHH 9X,3HL/0 7X,6HL(G)/L 5*,6HL('U!/- 

1 4X , 8HL AM BOA T 4X,5H?S/PA 8X,Z3VC 6 X , 3HGRF . C3EF AX, 6H1*D**2 5X» 

2 7HK = CP/C V ) 

IF (.NUT. VARA) WRITE (5,bl) 8ETAD, ALF A ,H , CL , LG , LU, LT , PS , VC , AO, D 1 , < 

6 1 FORMA T ( 12GL 1- 3) 

IF (VARA) WRITE (6,o2) BET AO , ALF A ,i) , CL , LG , LO, L T , PS ,VC , 01, < 

62 FORMAT ( 9G 1 1 . 0 , 8H VA R I A 8 LE OX 2311.3) 

ZG = LG/2. 

ZO = LO/2 . 

C2A = 2. 

IC3A = (0 .,2. )*ASC*H32/D2*L 
IC4A = IC3A 

C5A = (H 3 + AA*CJ*CO-M32) /D2*L*L 

IC6A = (0.,3.)#ASC*H31*(ALFA*(H2 — (-l.)-ALFI^H2*(H2+H-l.) )/ D2/DM*L 
C6B = L/C l 

C4B = — C6 8* ( ALF 1 + ALPA*H) 

CS5 = L*L 

D 3 T = 3./DN/02*ALFA*32*H-l.)*(l.+ALFI*(H3-l.)<‘SI*SI’M2.*ALFA 
1 + ALF I*(H3+1. ) ) ) 

DG = IC3A*Cl*L/2. 

DGP = C 1* L 

00 300 1=1, IS 

SOS = SIGN ( 1. , SL AM ( I ) ) 

SL AM( I ) = ABS ( SLAM! 1 ) ) 

LAMBDA = SL AM ( 1 ) 

IF (L AMBDA.EO.O..ANO. (PS.cw. 1 .. OR. LT.EU.O. ) ) GO TO 300 
IF (LAMBDA. EO.O.. AND. ( PS. E 0 . 1 . . Ok. l T. t Q. J . ) ) GO TO 300 
IF (LG .EO .0. .AND. I . GT » 1) 3U TU 27 
LS = SGN*L*LAMBDA 
CP = ASC*H3i*L S/0 2 
WRITE (6,8) LAMBDA 

8 FORMAT ( 17HLS3LUTI UN FOR PUC 8X 7HL AM6DA= G10.3/J 

HL M = 0. 

HLM 2 = 0. 

MO = 0. 

PC = 1 . + CP*ZG 
IF (LT.EQ.O.) GO TU 25 
PC = AMAX l (PC, SURT( PS) I 
XF = ZG 
COX = LG/32. 

DO 20 J = 1, 10 
IF (.NOT. VARA) GO TO 15 
P = PC /PS 

IF (P .GT. 1. ) P = 1. /P 
AU = ALPHA ( P ) 

15 MO = FMU ( PC , P S ,OMO) * AO 
OMO = AO* DMO 
OMP = DMO /PS 
UMO = 2,*Cl*CP/LT 



HLM = LT*AMIN1< MU, UMO) /2 . 

HLM 2 = HLM*L 
HlM = HLM /CL 
I F (MO.GE.UMU) GO TO 81 
IF (CP .£ J .0. ) GO TO 22 

POG = SQR T(AMAXl(PC*PC +• L* L T *MO* ( ZO-ZG ) ,0.)) 

PL = ALOG ( (P03*CP~HLM) /(CP-HLM) ) 

ZG1 = (PJG-L. +■ HLM/C P*PL ) /CP 
POG = POG*CP 

OZG = (PG + CL*LI*(Z0-Z3) *DMP) / ( POG-HLM) *■ L T / ( 2 . *C l*L P*CP ) *DM P 
L * l P L + HLM* (POG-GP) /( POG-HLM) / (CP-HLM) ) 

PC L = PC - (ZGl-ZGI/JZG 
GO TO 24 

22 IF (ZG.EO.O.) GO TO 28 

PG2 = PC*PC * L*LT*MO*( ZO-ZG) 

PCI = PC - (HLM*ZG-.5*.5*PG2) / ( PC* . 5 * ( l . - PC* PC ) / MO*DMP ) 

GO TO 24 

28 ZG1 = ZU - ( l.-PC*PC) /HLM2/2. 

OZG = PC /HLM 2 + ( ZO-ZGL) /MU*OMP 
PCI = PC - ZG 1 /OZG 

IF (PC .LT .PGA .AND. PCI .LT.PCA) PC 1 = SORT ( l. . *HL M2*Z 0) 

GO TO 24 
8 L VI 2 ) = i. 

IF (LG.EG.O.) GO TO 92 
Y ( l ) = 0. 

CALL RKGC (F1,Y,XF,CJX,A,3,2) 

92 PCI = SQRT(AMAX1(REAL(YI2) ) **2 ~ L *L T* AMI Nl < M0, J MO } *( Z 0-7. 3 ) , 0. ) ) 

24 WRITE (6,10) J , PG , PC 1 , MO , UMO , C , AO 

10 FORMAT (18, 4X 2HPC 315.7.3HPC1 G 15 . 7 ,2 HM0G12 . 4 , 5HM AX MOGl 2 . 4, A4 , A6 , 
1 3X 9 HOT IF GOEF 310.3) 

IF ( A 3 S( P C 1 /P C - L. i.Lfc.EPS) GO TO 25 
IF (PC .LT .PGA. ANJ. PCI .LT.PCA) PC=PCl 
IF (MJ.LT .UMO. AND. J.GT.l) PC = PCI 
20 PC = (PC + PCD/2. 

SWI I ) =0. 

GO TO 300 

25 PSO = -LI /2.*L)M0/PC 
PSl= VC*CL /PC 

PCS = PC /PS 
U=L*MO*LT/PC 
OD = 0 * ( D l + . 5 ) /L) l 
0 = U/PC/2. 

WRITE (6,26) PC ,PCS,PS0,PSl ,RKEP 

26 FORMAT ( 6HK PC G1 0. 3 , 58 PC/P S G10.3.6H PSO G10.3,6H PS! 

1 ' G 10 . 3 , 32H RKOX REOUCEO UNTIL ERROR .LE. GL0.3) 

PSO = L*P SO/PS 
PSl = L*P S l 

27 SPC( I ) = PC 
SMO ( I ) = MO 
SAO ( I ) = AO 
SIG1 = SS I G ( I ) 

OS = 1. 

IF (LAMBDA. NE.O.) OS = . 4* LAMBDA 
C2B = -CP 

C4C = AS*H31*L S*L/D2 

C6C = 3.* AS* (H- 1. ) **2*H2*LS/D2/DN*L 

C7A = 3.*ASC* (H-l. ) * l H2 - 1 . ) *L S* l H2 * (H2 + L. ) + AS *H3 2 ) / 02/ 0 2 
C8B = CP 
UP = DPT* HLM 2 
DC = C 7A* C 1* L 
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WRITE (6, 1*0) 

140 FORMAT ( 8HK. LAMdM 4X , 83 F R/E P ALO 3 X , 8HF T / E PAL D 4X, IH/if EPA_ I) 2< , 

1 4X 3HPHI 4X 9HMCW2/PALJ 7X 4HRKDX 3X 12HMPC57 2LR5MJ2 IX 

2 11HWHIKL RATIO IX 12HVIB NR SI G HA 2 X 8HPS0 COEF ) 

FIRSI = .TRUE. 

CALL HYST (O..GIJ) 

FIRST = .FALSE. 

IF (LAMBDA .EQ. 0. ) 310=.2S 

SW( I ) = w 

SPHU I ) =PHI 

CALL HYST (SIGl, Gil) 

SIG = SIG 1 

IF (AdS(3ll/3IO).LE.EPSw) GO TO 294 

051 = SDS* SIGN US, 311) 

SIG2 = SI G l 

DO ISO J = 1, 20 
SIG2 = SIG2 + D SI 
CALL HYST (SIG2, G12) 

SIG = SI32 

IF ( Ad S(G I 2/3 10) .LE.EPSW) 30 TO 299 

052 = SOS* SIGM US, 312) 

IF (DS2 - DS1) 220, 210, 220 
210 SIGl = SIG2 
150 Gil = GI2 
GO TO 299 

220 DO 235 J = l, 10 

SIG = (G12*SIGI - 311* SI 32) /(GI2 - G 1 1 ) 

SD(I) = ( G 12 - G11J/ISIG2 - SIGl) 

CALL HYST (SIG, GI) 

IF(ABS(GI/GIG).LE.EPSW) 30 TO 299 
IF (ABS(3I2).GT.4BS(3I 1) ) GO TO 230 
Gil = G I 2 
SIGl = SI G 2 
230 G 12 = GI 
235 SIG2 = SIG 

299 Sw 3( I ) =WR 
ACM ( I ) = CMA 
SCM ( I ) =CMC 

SS IGl I ) = SIG 
SPS( 1 ) = CPS 

300 COMTINUfc 
WRITE (6, 707) 

IF (SGN.LT.O.) WRITE (6,909) 

IF ( SGN.EQ.O. ) WRITE (6,390) 

IF (SGN.LT.O.) WR I TE (6,608) 

WRITE (6,600) CL, LO,L T,0i ,K 

600 FORMAT ( IHL 8X,3HL/D 8X,4HL0/L b X, 8 HL AMBOA T 4 X , 10 H 1+DELT A **2 3X, 

1 7HK=CP/CV 4X , 9HGR1 F COcF 2 X , l OHREC . VOL . VC /5X.5G12.3) 

IF (..MOT. VARA) WRITE (6,601) AO, VC 

601 FORMAT (1H+ 64X 2G12.3) 

IF (VARA) WRITE (6,602) VC 

602 FORMAT ( 1H+ 65X dH VARIABLE G15.3) 

WRITE (6, 603) P S, BE TAD, ALFA ,H ,LG ,RKE P 

603 FORMA T ( /8X 5HPS/PA 7X 4H8ETA 8X 4HALFA LIX IHH 10X 4HLG/L 8X 
1 4HRKEP /5X 7G12.3) 

WRITE (6, 32 5) { SL AM ( I ) , SW (I ) , SPH l ( D , SW3 ( I ) , S S I G( I ) ,ACM( I ) , S CM ( I ) , 
l SPC(I),SD(l), SMJ ( I ) , SAD ( I ) , S P'S ( I ) , I = 1, IS) 
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3 2b FORMAT (/1HL 3X.6HBKS NR 3X,8HRR3 LUAD 2X ,2 ZHATT ANGLE */HiR_ RAT 
LID 2X t faHY 1 8 NR 3X , 2 ON S T A 3 I L I I Y PARAMETERS liX 3HDIF1) 6 X '+HM AS S 
2 4X 7HORIFICE 3X 3HPSO / 

1 4X 6 HLAMBDA 4X 7HW/EPAlD 5X 3HPHI 6 X 8 HWP/UMEGA 4X 5HSIGMA 4X 

2 30 HMC W2/PALD MPC3/2LR3MU2 PC/PA 4X BrIL) (SIGMA) 3X 7HFLUW MO 

3 2X 7HC0EF AO 3X 4ri3iibF / 

4 ( / 3X 2G LO. 3jGll» '3 »t>12. 3» 2G 10. 3 9 20 L2»3» OR. 2»jL 1. 3 tOLOo '3 *iS. 3)) 

GO TO 3 

END 


ilBFTC HYDST 

SUBROUTINE HYST (SIGMA, 31) 

COMMON /C F F / HLMjCP ,G2A,C2B,I C3A , I C 4 A » 1C4W , ICS4W, IC6A, l Cs*l , ICSSRi 
L C5A» C 7A , C 83 , HLM2 ,C S5 

COMMON /CHS/ LAMBDA ,PSU ,PSl »C MA ,CMC , W, PHI » WR, RKFP, W, L, PC, C4B, C4C, 
1 C68.C6C ,FIR ST ,GI Z,3PS 

COMMON /CHR/ CG,OD,Z3,ZO , 

REAL LAMBDA,L 

COMPLEX AAjAB ,BA,BB,CA,CB,CG, TD X 

COMPLEX lC3A»IC4A»IC4WfIGS4W,IC6AslCbyj,ICSbW 

LOGICAL F IRST 

IF (F IRS I .OR . SIGMA. NE. 0. ) 30 TO 10 
Gl = GIZ 
WR = 1000 » 

IF (LAMBDA .Nfc . 0. ) WR = 0. 

CMC = 0. 

CM A = 0, 

RETURN 

10 CPS = l. 

CG = PSO*CPS +■ Q + (0. ,I.)*SIGMA*PSL 
C P SU = 0 MEANS ORIFICE MASS FLOW DOES NOT Y ARY WITH THCTA 
DX = .23 
TDX = 0. 

WR=L AM -3D A - SIGMA 

IC4W = (0.,L.)*(u4B*WR + C4C ) 

IC6W = ( D., V. }*(C68*WR + C6C) 

ICS6W = ( 0., L. ) *C S54WR 
ICS4W = - ICSfaW 
DO 110 N = l, B 
BA = 0 . 

BB = { -1. ,-l. » 

CALL RKLP { B A » DX , AA , 3 A ) 

CALL RKLP (BB»QX»AB,3B) 

CA = - 3. I 41392 63 * ( AB*CA - AA*CB)/(AB - AA ) 

IF (CABS( TDX/CA-1 . ) .Lt. RREP) GO 10 120 
DX = „3*DX 


36 



110 TOX = CA 
120 RDX=OX 

GR = REAL (CA) 

GI = A IMA GCA ) 

IF IF IRST ) GI Z = 31 
W = CABS(CA) 

PHI = ATAN2(GI ,GR )*S7. 2957796 
wR = LOOO. - 
CMC = 0. 

IF (SIGMA. NE.O.) CMC = I 4 A. *GR/ S IGM A**2 
CMA = CMC *L AM BO A* *2 /3t». 

IF ( L AMBD A .Mb . 0. ) WR = . 5* 31 G MA / LAMB J A 

WRITfc (6,143) LAMBOA.GR, GI ,W, PHl,CM4»Rl)X»CMC»WR»SlGMA,CPS 
143 FORMAT (3G11.3, 3 i 2. 4 , 3 L 1 . 3 ,G 12. 4 , 2GI 3 . 4 , 3G1 2 . 4 ) 

RETURN 

END 


$ IBFTC FM 

FUNCTIUN FMO (PC, PS, OMO) 

COMMUN /C FM / PCR , K , SKI , SK2 , VARA , AU 
COMMON /C HK /CC , UC , C ( 2 ) 

REAL K 

LOGICAL VARA 

OALPH m = . 06431 ♦- (-2. *.76193 * 3. * .434 15 * V ) *Y 
IF {PC.GT.PS) GO TO 20 
PCS = PC/PS 
AU = OAL PHI PC S) 

PM = -PS 
GO TU 30 
20 PCS= PS/PC 

AO = -PCS*PCS*DALPH (PC S) 

PM = PC 

30 IF (PCS. GT. PCR) 30 TO 40 

FMO = PM* SKI 
DMO = 0. 

IF (PC.GT.PS) UMO = P 3* SK 1 
IF (VARA) UMO = OMO *■ F MO*AD / AO 
C( 1 ) = CC 
RETURN 

40 PK = PCS**! i. /K) 

PSJ = SORT (I.-PCS/PK) 

FMO = PM* SK2* PK*P SO 

IF (PS.GT.PC) OMO = -PS*SK2*( PK/PCS- (K+l. )/2. >/ PSQ/K 

IF (PS.LT.PC) OMO = PS*SK2*(K-l.) / K* { PK-PCS /2 . ) / PS 0 

IF (PS.EJ.PC) OMO = 1.E10 

IF (VARA) OMO = OMO * FMU*AU/AO 

C( 1 )=UC 

RETURN 

ENO 
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$IBFTC BU 


BLOCK. DATA 

COMMON /OHK/CC » UC *0(2 ) 

DATA CC,UC»C<2) / 1H ,40 UN • 6HCH0KE0 / 
END 


$ I BFTC RKLPS 

SUBROUTINE RKLP ( G I N , DX , AD , CD ) 

COMMON /CY4/ UP , DC , OG , JG P 
COMMON /CHR/ CG,0D,Z3,ZU 

COMPLEX Y ( 5 ) , A { t» J ,B(5) ,G I N , AD ,CI3 ,C3 , XF,COX 
COMPLEX DG 
EXTERNAL F2,P3,F4 
Y( 1 1= 0. 

Y( 2) = I. 

Y( 31 « 0. 

Y< 4) = GIN 
Y< 5) =0. 

IF (ZG.EU.O.) GO TO 10 
XF = ZG 

COX = ZG/AMAXO(lFIX(ZG/DXf. 1) ,11 
CAlL RKGC <F2, Y,XF ,COX,A ,B ,5 ) 

Y ( 4 ) = DP/Y12I + DC + UG<= Y (3 ) * DGP*Y(4) 

10 IF (ZG.EQ.ZG) GO TO LB 

XF = ZO 

CDX * (ZU-ZG) /AMAXOI I FI X( (ZO-ZG) /DX+. I ) ,1 ) 
CALL RKGC (F3, Y,XF ,CJX,A,3,5 I 
15 Y( 4) = Y( 4) - QD - C3*Y(3) 

IF (ZO .EQ..5) GO TO 20 
XF = .5 

CDX = (.5-Z0)/AMAX0(IFIX( (.5-Z0J/DX+.I) ,1) 
CALL RKGC ( F 4 , Y , XF ,0 0 X , A , B , 5 ) 

20 AD = Y ( 4) 

CD * Y { 5) 

RETURN 

END 
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SIBFTC FF 


SUBROUTINE FI (Y»YP) 

COMMON /C FF/ HLM.CP ,C2A,C2b, IC3A.LC44, IC4W, LCS4W, IC6A, IC6X , IC36W, 
1 C5A,C 7A,C 33 ,HLM2 ,C S5 

COMPLE X IC3A , IC4A » !C4w,lC S4 W, I C6 A, I C6 w , ICS5R f Y ( 1) , YP{ 1 ) 

C SUBROUTINE FOR PU 
YP I 1) = 1. 

YPt 2) = CP- HLM/YI2) 

RETURN 

ENTRY F2 (Y,YP) 

C SUBROUTINE FOR G IN GROUVED REG I ON 
YP{ 2) = CP -HL T/YI2) 

YP ( 4) = -( <C2A*YP ( 2UC2B) /Y(2) + IC3A )*Y<4) - (t 

1 IC4A*YP ( 2H-IC4W *■ I-YP(2)+C8B)*YP(2)/Y{2) ) / Y ( 2) -C5A )*Y( 3) 

2 + IC6A*YP(2)+IC6W - C7A*YP(2) /Y<2) 

GO TO 10 

ENTRY F 3 (Y,YP) 

C SUBROUTINE FOR G IN SMOG TO REGION 
YP( 2) = -HLM2/YI2) 

YP I 4) = - C 2A* YP ( 2 ) * YI 4 ) / Y (2) - ( ( I CS4W - YP ( 2 ) *YP( 2 ) 

1 /YI2I1/YI2) - C$5 ) *Y I 3) *- 1CS6W 
GU TO 10 
ENTRY F4 < Y, YP ) 

C SUBROUTINE FOR G IN REGION INBQARU OF ORIFICES 
YP( 2) = 0. 

YP ( 4) = (CS5 - 1CS4W/YI2J ) # Y( 3 ) +• ICSoW 
10 YP ( 1 ) = 1. 

YP ( 3) = V ( 4) 

YPt 5) = Y( 3) 

RETURN 

END 


i IBFTC RKGHC 

SUBROUTINE RKGC (OERIV, Y, XFINAL, DELTA, U, YP, N ) 

C RUNGE-KUTTA-GIlL INTEGRATION OF N-l COMPLEX FUNCTIONS Of COMPLEX 
C ARGUMENT Y( 1 ) FROM Ytl) INI UAL TU XFINAL 
DIMENSION A ( 4 ) , 3(4), C(4) 

COMPLEX Y ( l ) , Q ( 1 ) , YP ( 1 ) , XFINAL, DELTA, T 

DATA A/. 5, « 2928932 10, 1.70710678, .166666657/, B/2„, l., 1., 2./, 
1 C/.5 , .292673219, 1.70710678, .5/ 

KK = CABS( ( XF INAL-YI 1 )) /DELTA) + .01 
IF (KK.LE.O) GU TU 20 
DO 5 I = 1, N 
6 Q(1)=0. 

DO 10 K = 1, KK 
DO 10 J =1, 4 
CALL OERlVtY.YPJ 
DO LO 1 = 1, N 

T=A! J )*{ YPID-B { J)*Q( I ) ) 

Y ( I ) = Y(I)+DELTA*T 
LO Q( I )=0(I)+3.< ! T-C(J) , XYP(I ) 

20 RETURN 

END 
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